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Abstract A bright prominence associated with a coronal mass ejection (CME) 
was seen erupting from the Sun on 9 April 2008. This prominence was tracked by 
both the Solar Terrestrial Relations Observatory (STEREO) EUVI and COR1 
telescopes, and was seen to rotate about the line of sight as it erupted; therefore, 
the event has been nicknamed the “Cartwheel CME.” The threads of the promi- 
nence in the core of the CME quite clearly indicate the structure of a weakly to 
moderately twisted flux rope throughout the field of view, up to heliocentric 
heights of 4 solar radii. Although the STEREO separation was 48°, it was 
possible to match some sharp features in the later part of the eruption as seen 
in the 304 A line in EUVI and in the Ha-sensitive bandpass of COR1 by both 
STEREO Ahead and Behind. These features could then be traced out in three- 
dimensional space, and reprojected into a view in which the eruption is directed 
towards the observer. The reconstructed view shows that the alignment of the 
prominence to the vertical axis rotates as it rises up to a leading-edge height of 
~2.5 solar radii, and then remains approximately constant. The alignment at 2.5 
solar radii differs by about 115° from the original filament orientation inferred 
from Ha and EUV data, and the height profile of the rotation, obtained here 
for the first time, shows that two thirds of the total rotation is reached within 
« 0.5 solar radii above the photosphere. These features are well reproduced by 
numerical simulations of an unstable moderately twisted flux rope embedded in 
external flux with a relatively strong shear field component. 
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1. Introduction 

The bandpass of the inner coronagraph (COR1) telescopes on the two Solar 
Terrestrial Relations Observatory (STEREO) spacecraft runs from 650-670 nm 
(Howard et al. , 2008). This range was selected to include the hydrogen Ha line 
at 656 nm, so as to be sensitive to erupting prominences associated with coronal 
mass ejections (CMEs). A particularly bright prominence eruption was observed 
by both COR1 telescopes on 9 April 2008. Although the relative contributions 
of Ha emissions and Thomson-scattered light cannot be derived from the COR1 
observations alone, it is presumed that Ha is a significant contributor for the 
extremely bright parts of the prominence early in the eruption. This CME and 
prominence eruption was also observed by the STEREO Extreme Ultraviolet 
Imager (EUVI) telescopes (Howard et al. , 2008) before entering the COR1 field 
of view, as well as by TRACE at 171 A (Handy et al ., 1999), and by the Hinode 
XRT telescope (Golub et a/., 2007) and EIS spectrometer (Culhane et al ., 2007). 
Landi et al. (2010) analyze this event, combining Hinode /XRT and EIS data, 
together with Solar and Heliospheric Observatory (SOHO)/UVCS/EIT/LASCO 
and STEREO/EUVI/COR1/COR2 images, to characterize the thermal proper- 
ties of the ejected plasma, and constrain the heating rate. The post-CME current 
sheet following this event is analyzed by Savage et al. (2010) and Ko et al. (2010). 
Patsourakos and Vourlidas (2011) determine the geometrical parameters of the 
CME and post-CME current sheet in the field of view of the outer coronagraph 
(COR2) on STEREO. This event has become known as the “Cartwheel CME” 
because of the highly visible rotation around the line of sight seen during the 
initial stages of the eruption. Rotation about the vertical direction can also be 
observed through triangulation as the prominence material rises through the 
EUVI and COR1 height ranges. The plane of sky motions which give this event 
its “Cartwheel” name is actually a combination of the rotation about the vertical 
axis together with a deflection of the prominence in latitude and longitude during 
the initial rise phase. A rotation around the axis of the prominence, known as 
“roll effect” (Martin, 2003), may have contributed as well. Further rotation about 
the vertical direction is indicated by the CME orientation in the COR2 data. 
This paper presents an exploration of the rotation about the vertical direction in 
the EUVI and COR1 height ranges, using triangulation from the two STEREO 
viewpoints. The height profile of the rotation angle for the erupting prominence 
is obtained, which could not be done prior to the STEREO mission. 

The separation of the two STEREO spacecraft at this time was 48°. Co- 
identification of features in both views is quite difficult with such a large sep- 
aration. However, it was possible to match some sharp features in the later 
part of the eruption as seen by EUVI in the He II line at 304 A. Features were 
also tracked in white light (with a possible Ha component) as seen by COR1, 
although not necessarily the same features seen at 304 A. 

The derived rotation and rise profiles of the prominence, as well as the 
STEREO images, are compared with the corresponding data obtained in a series 
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Figure 1. Ha image of the filament as seen on the disk on 5 April 2008, 01:51 UT, before the 
first eruption on that day, and on 6 April, 02:44 UT (circled), reformed after the two eruptions 
on 5 April (from the Yunnan Astronomical Observatory in Kunming, China). 


of CME simulations in a companion paper (Kliem, Torok, and Thompson, 2011; 
in the following Paper II), in order to find the causes of the large total rotation 
observed. The comparison suggests that the major part of the rotation was due 
to the presence of a shear field component of considerable strength in the source 
region (Isenberg and Forbes, 2007) and that the weak helical kink instability of 
a moderately twisted flux rope (Torok, Kliem, and Titov, 2004) contributed as 
well. Field line plots of the best matching case are included here to demonstrate 
the correspondence with the observed overall shape and helical threads of the 
rising prominence. 


2. Data Analysis 

2.1. Source Region 

The eruption on 9 April 2008 occurred in the remnants of NO A A active region 
(AR) 10989, located about 23° behind the west limb as seen from Earth, and 
had an onset time near 09 UT. The filament had been visible nearly through- 
out the disk passage of the active region, which was in its decaying stage and 
spotless after 31 March. The filament erupted twice on 5 April 2008, reforming 
afterwards. These eruptions were visable on disk to EUVI on STEREO Ahead. 
Both eruptions appeared to propagate inclined from the radial direction toward 
higher southern latitudes, similar to the eruption on 9 April. The inclination 
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Figure 2. STEREO Ahead EUVI 195 A images of AR 10989 before and after the second 
filament eruption on 5 April 2008, displaying an area of 640" on a side. The image on the left 
at 19:45:30 UT shows the full extent of the filament on the souteastern end at the onset of 
the eruption (up to the middle of the circle), and the image on the right at 21:05:30 UT (with 
the circle plotted at the same position) shows that a dimming has developed at this location. 
Therefore, the erupting flux should be rooted in this area. Comparison with Fig. 1 shows that 
only the upper branch of the filament is visible in Ha. 



Figure 3. STEREO Ahead EUVI 171 A image of the reformed filament on 6 April 2008, 
15:01 UT. The area shown is 320" on a side. The circles mark the ends of the filament material 
visible at this wavelength. Their distance is « 175 Mm. The main body of the filament as seen 
in this line is oriented at an inclination of £s26° to the east- west direction on the Sun. 




Figure 4. Magnetogram of AR 10989 on 31 March (left) and on 4 April 2008 (right) taken 
by SOHO/MDI. The overplotted bars indicate rough estimates of the distance between the 
center of each polarity (in a center-of-gravity sense), which is an important input parameter 
for the numerical modeling of the event in Paper II. 
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made it difficult to discern rotational motions based on images from only a 
single viewpoint. However, a stereoscopic reconstruction of the second eruption 
revealed a considerable rotation of about 90° (Bi et al ., 2011). Interestingly, 
the rotation was in the clockwise direction, opposite to that of the eruption on 
9 April. 

As with most filaments, it is difficult to discern the magnetic connections 
at its ends. However, both eruptions on 5 April produced a dimming and end- 
point brightenings in the EUV (Wang, Muglach, and Kliem, 2009) near the 
southeastern end of the filament, which were located in negative polarity, the 
polarity which dominated the southwestern side of the filament channel, so that 
the filament was dextral according to the classification of Martin (1998). This 
is further supported by the right-bearing orientation of the filament barbs (see 
Figure 1). Dextral filaments are embedded in field of left handed chirality, as 
indicated by the typical skew of the overlying arcade of coronal loops (Martin, 
1998). If the assumption of flux rope topology holds for the filament, then the 
filament itself is also threaded by left-handed field. Otherwise, it may have the 
opposite chirality (see Martin and McAllister, 1997; Ruzmaikin, Martin, and Hu, 
2003; Muglach, Wang, and Kliem, 2009). If a filament in dominantly left-handed 
helical field rotates upon eruption, the direction of rotation is generally found 
to be in the counter-clockwise direction, as expected from helicity conservation 
(Rust and LaBonte, 2005; Green et a/., 2007). 

On 5 April, the filament was oriented at an inclination of ~ 24° to the east- 
west direction on the Sun and had a length of ~ 135 Mm in the image shown 
in Figure 1 (top). However, the dimmings and endpoint brightenings formed in 
the eruptions later on the same day suggest that the true extent of the flux in 
the filament channel was larger by about one third, ^175 Mm, see Figure 2. 

Due to foreshortening, the Ha images on the following day show the filament, 
its magnetic connections and barbs far less clearly, but it can be seen that the 
filament reformed in a similar location (Figure 1, bottom). The better perspec- 
tive of the EUVI- Ahead images indicates that the main body of the reformed 
filament was nearly straight, oriented at a tilt angle of « 26° to the east- west 
direction, and extended across a similar length as on 5 April (Figure 3). This is 
most clear if the 171 A and 304 A images are viewed in animated format. As 
far as the increasing foreshortening allows a judgment, the filament appeared to 
keep its orientation and shape through the subsequent days until a new eruption 
launched the Cartwheel CME. 

Since the numerical modeling in Paper II reveals a strong influence of the 
distance between the polarities in the source region on the rotation of the ejected 
flux rope, we attempt to estimate this parameter. Figure 4 shows magnetograms 
of the region from the SOHO Michelson Doppler Imager (MDI) on 31 March and 
4 April 2008. The overplotted bars run through the middle of each polarity in a 
one-dimensional center- of- gravity sense, providing an indication of the distance 
between the main polarities. This value increases from ^ 40 Mm to a range of 
^ (40-75) Mm over the four-day time span. A range twice as wide may have been 
characteristic of the configuration after the similar time span to the eruption on 
9 April. It appears impossible to estimate a single relevant value. We will adopt a 
distance of 90 Mm as a base value in most simulations in Paper II, but also study 
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Figure 5. Location of the heliospheric current sheet for Carrington Rotation 2068 from a 
potential-field source-surface extrapolation. The current sheet runs nearly east-west above the 
estimated CME direction of propagation as indicated by the vertical bar. The green and red 
regions represent open field regions of negative and positive polarity respectively. 


the influence of its variation. This base value lies in the middle of the estimated 
range, so it is consistent with the fact that much of the rotation occurred in 
the course of the radial ascent above the south-east edge of the remnant active 
region, which is characterized by a considerable distance between the polarities 
(see Section 2.3). 

We have also estimated the height of the prominence from data taken by 
the EUV Imaging Telescope (EIT) (Delaboudiniere et al . , 1995) onboard SOHO 
on 8 April, when the prominence was close to the limb. The material seen at 
that time in the 304 A line extended up to a height of « 0.04 R © above the 
photosphere, and the absorbing material seen at 171 and 195 A extended up 
to ~ 0.033 Rq. From EUVI 304 A data about an hour before the eruption, we 
estimate the height to lie in the range « 0.05-0.06 i^©, so it is clear that the 
prominence experienced a slow rise before the eruption. 

Motions along the prominence could be seen in the 171 A channel of EUVI- 
Ahead on 9 April from about 08:20 UT onward. The first upward motions of 
prominence material along the CME path can be discerned by comparing the 
images at 08:48:30 and 08:51 UT. Landi et al (2010) give a rather conservative 
estimate of the CME start time, 09:10 UT, but the 171 A images, taken at 2.5 min 
cadence, show that most of the prominence body visible to EUVI-Ahead was 
already moving by 08:53:30 UT, the start time quoted in Savage et al. (2010). 
Since the event commenced in the absence of a strong perturbation (no signs 
of a significant brightening, of a jet, or of perturbations resulting from nearby 
activity were seen), it must have developed from a small perturbation when the 
configuration was near the boundary between stable and unstable states. In such 
cases the initial motion of the unstable flux is expected to behave exponentially, 
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which requires some time to develop to a level that causes visible changes. Rapid 
changes in the vertical position were first seen between the EUVI- Ahead images 
at 08:48:30 and 08:51 UT. Therefore, the actual start time should lie before 
08:51 UT. 

A helmet streamer, best seen in the COR2-Ahead images, extended above 
the active region. In order to estimate the orientation of the heliospheric current 
sheet above the streamer, we ran a potential- field source-surface extrapolation 
(PFSS) and a Wang-Sheeley-Arge model for Carrington Rotation 2068 at the 
Community Coordinated Modeling Center (CCMC), using a source surface lo- 
cation at 2.5 solar radii for both models. Figure 5 shows the output of the PFSS 
model, which yields a nearly east-west orientation of the heliospheric current 
sheet above AR 10989. The Wang-Sheeley-Arge model gives the same result. 

2.2. 3D Reconstruction 

The SolarSoft routine scc_measure was used to measure the three-dimensional 
location of features in the STEREO EUVI and COR1 images. The action of 
this program is as follows: The user is presented with two side-by-side images, 
one from each of the two STEREO spacecraft. The images are selected so that 
they represent the same observation time, with a slight offset to account for 
the difference in light travel time from the Sun. The user can zoom in on the 
region of interest in the two images, and adjust the color table and data range 
to optimize the appearance of the feature being measured. A point is selected on 
one image with the cursor. The program calculates the three-dimensional line 
of sight represented by this point, and then overplots the projection of this line 
onto the image from the other satellite. This is known as an epipolar line. Since 
both the EUVI and COR1 optics produce a gnomonic projection on the CCD 
detector, straight lines in space will always appear as straight lines in the image. 
The feature selected by the user in the first image must appear along the epipolar 
line drawn by the program in the second image. The correct location along this 
line is selected by the user, which leads to another line of sight calculation which 
intersects the original line of sight. The intersection of these two lines determines 
the three-dimensional (3D) location of the feature. 

The main difficulty in applying this technique to the prominence eruption is 
source confusion. On 9 April 2008, the two STEREO spacecraft were separated 
by 48° along the ecliptic plane. Thus, the appearance seen by STEREO Ahead 
was quite different from that seen by STEREO Behind. This made it very difficult 
to locate features which could be positively identified to be the same in both 
images. Another goal was to identify features which could be tracked through 
several frames. (However, no attempt was made to match identifications between 
EUVI and COR1.) In spite of these difficulties, it was possible to identify several 
features which could be identified in both views, and which could be tracked 
through several frames. 

2.3. EUVI and COR Measurements 

Figures 6 and 7 show the tracked features as seen in the 304 A channel of EUVI- 
Ahead. The He II 304 A channel was chosen as being most representative of 
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Figure 6. Images and plots of the prominence eruption from 10:06 to 10:26 UT, as seen in 
the 304 A channel of the STEREO EUVI telescopes. Subsequent frames are shown in Figure 7. 
Along the left are shown the EUVI-Ahead images with the tracked features overplotted. The 
plots along the right show the same data reprojected to a viewpoint at Stonyhurst longitude 
98° west (relative to Earth) and latitude 24° south. Axes are in units of solar radii. 


the cool prominence material seen later in a combination of Ha and Thomson- 
scattered light by COR1. Although the eruption is clearly seen by Ahead at 
9:06 UT, and by Behind at 9:36 UT, source confusion made it impossible to 
track features earlier than 10:06 UT. Initially, little sense could be made from 
an examination of these data viewed in 3D (e.g., through an anaglyph rep- 
resentation). However, because the measurements are made in 3D, they can be 
reprojected to another viewpoint. By selecting the proper viewpoint, it’s possible 
to see aspects of the data that are not evident in the frame in which the data 
are taken. The plots to the right in Figures 6 and 7 show the same data as in 
the images, but from a viewpoint at Stonyhurst longtitude 98° west (relative 
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Figure 7. Continuation of Figure 6 for time steps 10:36 to 10:56 UT. 


to Earth), and 24° south. It is estimated that from this viewpoint, the CME is 
traveling straight toward the observer. 

From this perspective, a certain organization can be perceived. The primary 
alignment in the 10:06 UT data is vertical, with the upper prominence threads 
lying almost directly over the lower threads. However, in the later images the 
primary alignment is distinctly sloped, with the upper threads lying toward 
the east, and the lower ones toward the west. This is particularly well seen in 
the image at 10:46 UT in Figure 7, where the data display an overall slope of 
—50°. As one goes through the time steps, a counter-clockwise rotation can be 
perceived as one steps from 10:06 to 10:46 UT. 

The prominence structure must have rotated through ~ 65° from its original 
orientation to achieve the nearly vertical structure seen at 10:06 UT, followed 
by an additional « 50° to reach the maximum rotation at 10:46 UT. Thus, the 
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full amount of rotation experienced by the prominence in the EUVI height range 
was on the order of 115°. 

The EUV observations of the early part of the eruption also show a shift of 
the entire prominence structure southward, from about -14° to -24° latitude. 
Compared with the earlier Ha observations (Figure 1), the prominence must 
have also shifted ^ 15° eastward in the early phase of the eruption. These 
initial motions were followed by radial propagation at the new latitude and 
longitude, ~98W24S. This is consistent with the analysis of Savage et al. (2010) 
and Patsourakos and Vourlidas (2011). The latter authors fitted a croissant- 
shaped flux rope model to COR2 images from both STEREO satellites and 
found that the radial propagation of the CME projects back to a position 11° ±5° 
behind the limb, as seen from Earth, and 17° db 3° south. Savage et al. (2010) 
quote a heliographic position of 113° for the associated active region at the 
time of the eruption, and describe the initial motion of the prominence as being 
toward the observer. The measurements shown in Figures 6 and 7 were derived 
from the radial phase of the eruption. 

The same analysis performed on EUVI was also performed on the COR1 
images. The results are shown in Figure 8. The prominence is seen to maintain 
the orientation seen in the final EUVI images. There is a slight indication of a 
shift back in the clockwise direction. However, the orientations for the final two 
time steps at 11:15 and 11:25 UT are based on only parts of the prominence 
legs, so the apparent reverse rotation may not truly represent the behavior of 
the loop as a whole. 

Patsourakos and Vourlidas (2011) estimated the orientation of the CME flux 
rope from their croissant model fit at a leading-edge heliocentric distance of 
about 13 Rq. They found that the flux rope (croissant) axis was inclined to 
the east- west direction on the Sun by a tilt angle of —4° db 7°. This indicates a 
further counter-clockwise rotation to a total value of 150° ± 7° from the original 
orientation estimated in Section 2.1. (The other possible interpretation of a 
clockwise rotation by « —145° following the counter-clockwise rotation in the 
EUVI height range appears far less likely.) It should be noted that the croissant 
model does not include any writhe of its axis, so that the uncertainties may be 
higher; however, this is not expected to change the result by a large amount. 
The rotation angle at large heights represents a close alignment of the CME 
flux rope with the heliospheric current sheet (Figure 5), which corresponds to 
the suggestion in Yurchyshyn (2008) and Yurchyshyn, Abramenko, and Tripathi 
(2009). 

To our knowledge, the rotation of the prominence by « 115° in the corona 
and by « 150° up to 13 Rq belongs to the largest ever inferred in this height 
range. Yurchyshyn, Abramenko, and Tripathi (2009) report one case of rotation 
by « 143° and six cases in the range 80°-100° out of a sample of 101 halo CMEs 
observed by SOHO/LASCO, z.e., up to a distance of 30 Rq. The largest rotation 
found at 1 AU is on the order of 160° (Dasso et al . , 2007; Harra et al . , 2007). 

It is also of interest to note that the COR1 images resolve the basic structure 
of the dense prominence material in the core of the CME out to a heliocentric 
distance of four solar radii. The CME core has the structure of a weakly to mod- 
erately twisted flux rope: a single flux loop with threads that are systematically 
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Figure 8. Images and plots of the prominence eruption from 10:45 to 11:25 UT, as seen in 
the COR1 telescopes. Along the left are shown the CORl-Ahead images with the tracked 
features overplotted. The plots along the right show the same data reprojected to a viewpoint 
at Stony hurst longitude 98° west (relative to Earth) and latitude 24° south. Axes are in units 
of solar radii. 
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Figure 9. Orientation of the prominence body as a function of time. The dashed line rep- 
resents the original orientation as estimated from the EUVI-Ahead 171 A image on 6 April 
in Figure 3. The points at 10:16 and 10:26 UT are considered to be less reliable than the 
others — see text. 


but only weakly to moderately inclined to the axis of the loop. These indications 
of twist prove to be the primary observational finding that sets a preference for 
a slightly kink-unstable flux rope above a kink- stable flux rope, which match 
the observed rise and rotation characteristics to a comparable degree in the 
numerical modeling of the event in Paper II. Only few observations have revealed 
a flux rope structure for an erupting prominence so clearly in this height range 
(see, e.g., Plunkett et a/., 2000). The images also indicate a nearly self-similar 
evolution of the flux rope throughout the instrument’s field of view of (1.5-4) R 0 . 
The outer coronagraph COR2 imaged the CME core as well but did not resolve 
the details of its structure. 

2.4. Rotation and Rise Profile 

Figure 9 summarizes the rotation of the erupting prominence as a function of 
time from both EUVI and COR1. The values plotted were obtained by aver- 
aging the orientations of the individual threads whose 3D positions could be 
reconstructed. The values obtained at 10:16 and 10:26 UT are considered to be 
less reliable than the others in the figure, since at these times the threads which 
could be located in 3D do not appear to be organized in a well defined structure, 
as, for example, flux rope legs (see the right panels in Figure 6). 

We also used scc_measure to derive a time-height curve for the prominence 
eruption. This was done by matching points at the leading edge of the promi- 
nence in both images. Since these points in the two images are not guaranteed 


SOLA: prominence_v3.1f.tex; 10 June 2011; 16:52; p. 12 


X EUVI 


n — 1 — 1 — 1 — i — 1 — 1 — 1 — i — 1 — 1 — 1 — r 


n — 1 — 1 — 1 — r 


~o 


o 3b +C0R1 


g> 

'CD 
JZ 

CD 
CD 
"O 

CD 2 

CD 

c 

T 5 
03 
CD 


- X 


X 


+- 


+ 


+ 


+ 


X 


A 

X 


09:40 10:00 10:20 10:40 11:00 11:20 11:40 

Time (UTC on 9-Apr-2008) 


Figure 10. Heliocentric heights of the prominence leading edge vs. time. The value for 
10:26 UT is shown as a lower limit, since the prominence is right at the edge of the EUVI- Ahead 
field-of-view at this time. 


to represent exactly the same point in 3D space, the results can only be treated 
as an approximation. However, since the extent of the prominence material is 
relatively small, the approximation is quite valid. The resulting time-height curve 
for the top of the prominence is shown in Figure 10. The COR1 values appear to 
be a bit higher than one would expect from an extrapolation of the EUVI data. 
That may represent a difference between the appearance of the prominence at 
304 A versus white light. In addition, as is obvious from Figure 6, the heights 
derived from EUVI data after 10:06 UT may systematically fall short of the true 
heights, since the traceable threads did not appear to extend up to the true 
top of the structure, which is not the case for the COR1 measurements. The 
estimated height at 10:26 UT may additionally be influenced by the edge of the 
EUVI- Ahead field-of-view. 

A quadratic fit to the EUVI data gives an initial speed of 67 km s _1 , accelerat- 
ing to 173 km s -1 as the prominence leaves the EUVI field-of-view. A quadratic 
fit to the heights measured by COR1 gives a velocity starting at 206 km s _1 , 
and accelerating to 379 km s -1 by the end of the COR1 data. There’s some 
evidence from the COR1 quadratic fit that the acceleration strongly decreased 
when the leading edge of the prominence reached a heliocentric height of 3 solar 
radii. This is borne out by a linear fit to the COR1 data points above 3 solar 
radii, which gives a slightly lower velocity of 327 zb 9.2 km s -1 . The combined 
EUVI and COR1 data in Figure 10 give the visual impression that much of the 
acceleration actually occurred up to « 2.5 solar radii. The COR2 data give a 
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Figure 11. Prominence rotation vs. heliocentric height of its leading edge. The diamond 
symbols represent EUVI and COR1 height data from Figure 10 interpolated to the times of 
EUVI orientation measurements in Figure 9. 


speed close to 400 km s _1 . All of these values refer to estimated positions of the 
apex point of the visible flux rope structure in the core of the CME. 

Landi et al. (2010) also derive a time-height profile for the leading edge of the 
prominence material, using a combination of STEREO Ahead and SOHO data. 
Their results are very close to those presented here, with the possible exception 
of the amount of acceleration. Landi et al. quote an acceleration of 59.6 m s -2 , 
while our results are more consistent with a smaller value of 37.9 db 4.0 m s -2 
over the time range covered by Figure 10. They also give the projected velocity 
of the CME leading edge in the plane of the sky of STEREO Ahead , which 
approaches 700 km s -1 , still rising, at a heliocentric height of 3 Rq. 

Combining the information in Figures 9 and 10, we obtain the rotation of the 
prominence as a function of its leading edge height, see Figure 11. (As there are 
only three time steps in common between the EUVI data sets, we’ve interpolated 
the EUVI and COR1 height data from Figure 10 to the times of the EUVI 
measurements in Figure 9 to better show the trend.) These data emphasize 
again that much of the rotation is acquired in the height range h < Rq/ 2, or 
initial footpoint separations, above the photosphere. 


3. Numerical Modeling 

Paper II presents a parametric study of various effects which lead to, or influence, 
the rotation of ejected flux ropes. The initial flux rope twist, the strength of the 
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Figure 12. Shape of the erupting flux in the best-matching numerical model — a weakly 
kink-unstable flux rope of initial twist 3.57 r — at various stages of the rise, compared with 
STEREO Ahead images at 10:16, 10:26, 10:55, and 11:25 UT. A set of field lines enclosing 
the magnetic axis of the simulated flux rope is plotted at the times the rope apex has reached 
the corresponding heights, h = 7.3, 9.6, 20, and 30ho (ho - initial apex height), which are 
obtained from the scaling of the footpoint distance, 3.3ho, to the value of 175 Mm estimated 
in Section 2.1. The observed and simulated flux rope are displayed at the same perspective: 
the line connecting the footpoints of the rope makes an angle of 26° with the line of sight 
and the vertical axis in the simulation is tilted away from the observer by 8° (so that the 
magnetogram, B z (x,y , 0, t) , is seen from the bottom side). The dotted line indicates the edge 
of the COR1 occulting disk. Since the visible threads in the EUVI images do not outline the 
complete shape of the flux rope, sections of the field lines in the corresponding simulation 
snapshots are highlighted for better visual correspondence. For each field line, 200 segments 
at the top and 150 segments at each bottom end are plotted with reduced line width. The 
segments are obtained from a numerical integration along the field line with adaptive step size 
and are all different, resulting in different lengths of the highlighted sections, analogous to the 
observations. 
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external shear field component, and the height profile of the overlying potential 
field are varied. The latter is related to the distance between the main polarities 
in the photosphere. Figure 12 shows field line plots from the numerical model 
that best fits the observations in their entirety (rise and rotation profiles and the 
STEREO images). In this model, a moderate initial twist of = 3.57 r triggers a 
weak helical kink instability, which lifts the initially force-free, toroidal flux rope 
into the torus-unstable range of heights and also contributes to the rotation. The 
ejection is primarily driven by the torus instability (Kliem and Torok, 2006), so 
that the model belongs to the loss-of-equilibrium category (Priest and Forbes, 
2002). The major part of the rotation in this run is due to the presence of 
an external shear field component (pointing along the prominence and polarity 
inversion line and due to sources outside the flux rope), as first suggested in 
Isenberg and Forbes (2007). The strength of the shear field at the initial flux 
rope apex is 2/3 of the external poloidal field which holds the rope in equilibrium. 

The overall shapes of the flux rope in the considerable height range included 
in Figure 12, as well as the inclination of the still weakly twisted field lines to the 
axis of the rope at the largest height, match the observations reasonably well. 
(Note that the rope has acquired a large part of its total rotation already at the 
first time selected; see Figure 9.) Similar agreement is demonstrated in Paper II 
for the height-rotation and time- height profiles. 

A model with subcritical initial twist, <f> = 2.57T, which allows only for the 
development of the torus instability, reproduces the rise profile similarly well, 
while the match with the rotation profile is somewhat worse. The field lines at 
the strongly expanded stage of the final COR1 image in Figure 12 are rather 
straight and hardly show any indication of the observed twist. Moreover, this 
flux rope requires a considerable initial perturbation to reach the torus-unstable 
range of heights, which is not supported by the pre-eruption height estimate and 
the initial dynamics of the prominence described in Section 2 (see Paper II). 


4. Summary and Conclusions 

The images of the STEREO EUVI and COR1 telescopes resolve the basic struc- 
ture of the erupting prominence in the core of the “Cartwheel CME” on 9 April 
2008 out to heliocentric distances of 4 Rq. A single flux rope with indications 
of weak to moderate twist is revealed. The flux rope expands approximately 
self-similar ly in the COR1 field of view of (1.5-4) R 0 . 

The true path of the prominence could be reconstructed from stereoscopic 
observations even when the angular separation between the STEREO spacecraft 
was as much as 48°. The reconstruction reveals that the prominence rotated 
counter-clockwise by a large angle of ^115° up to a heliocentric height of 2.5 R 0 , 
where the rotation leveled off. A gentle backward rotation by ~ 15° may have 
followed in the height range up to 3.3 Rq. Two thirds of the rotation were 
acquired within 0.5 Rq from the photosphere. The coronal height profile of the 
rotation angle in a CME is thus derived for the first time. To our knowledge, 
this ranks as one of the largest rotations so far measured in the corona. 

Taken jointly with recent results about the orientation of the CME at a helio- 
centric distance of ~ 13 R 0 , obtained from fitting a croissant-shaped flux rope 
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model to stereoscopic images from C0R2 (Patsourakos and Vourlidas, 2011), a 
further counter-clockwise rotation by ~35° from the value at 2.5 Rq is indicated. 
This aligned the erupted flux closely with the heliospheric current sheet above 
the streamer extending from the active region. 

The rotation and rise profiles of the prominence in the EUVI-COR1 height 
range are reproduced by a numerical model that follows the evolution of an 
unstable force-free flux rope. A parametric study of flux rope rotation in this 
model, detailed in the companion paper (Kliem, Tdrok, and Thompson, 2011), 
suggests that the main part of the rotation in the Cartwheel event was caused 
by the shear field in the source volume. While the rotation and rise profiles 
can be modeled nearly equally well by a weakly kink-unstable flux rope of 3.57T 
initial twist and by a kink-stable flux rope of 2.57T initial twist, the indications 
of twist in the COR1 images and the absence of a strong initial perturbation 
in the EUVI- Ahead data favor the 3.57T model, suggesting that the helical kink 
instability did occur and contributed the remaining part of the total rotation. 
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